Condensation Energy of a Spin-1/2 Strongly Interacting Fermi Gas 
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We report a measurement of the condensation energy of a two-component Fermi gas with tunable 
interactions. From the equation of state of the gas, we infer the properties of the normal phase in 
the zero-temperature limit. By comparing the pressure of the normal phase at T = to that of the 
low-temperature superfiuid phase, we deduce the condensation energy, i.e. the energy gain of the 
system in being in the superfiuid rather than normal state. We compare our measurements to a 
ladder approximation description of the normal phase, and to a fixed node Monte-Carlo approach, 
finding excellent agreement. We discuss the relationship between condensation energy and pairing 
gap in the BEC-BCS crossover. 
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From a thermodynamic point of view, a superconduct- 
ing state is favored compared to a normal state when the 
free energy of the former (Eg) is lower than the latter 
(En). This energy difference, called condensation energy 
is a central concept in the BCS theory of conventional su- 
perconductivity. For example, in the weakly interacting 
regime the condensation energy is related to the super- 
fluid pairing gap A by 
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where Nf is the density of states at the Fermi energy pQ . 
For superconductors, the condensation energy is obtained 
from the measurement of the critical magnetic field H c 
at which superconductivity is quenched 
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where no is the vacuum magnetic permeability [ 1 . . While 
BCS theory (and relation ([I])) have proven very success- 
ful to explain conventional superconductivity, a similar 
description to explain exotic forms of superconductivity, 
such as encountered in cuprate or iron-compound mate- 
rials is still lacking. In particular, the role of the conden- 
sation energy in high-T c superconductors is thought to 
give insight on the mechanism that could be responsible 
for driving the superconducting transition (see for exam- 
ple [5H5] and references therein), though the method to 
extract it from experimental data is still a hotly debated 
issue [MS]- 

Ultracold atoms are now increasingly used as testbeds 
to experimentally explore quantum many-body physics, 



'Present email and address: nn270@cam.ac.uk, Cavendish Labo- 
ratory, University of Cambridge, J.J. Thomson Ave., Cambridge 
CB3 OHE, United Kingdom 
tleyronas@lps.ens.fr 



owing to their high degree of control [9] . It now becomes 
possible to simulate hamiltonians from various fields of 
physics, such as neutron matter or condensed matter 
physics in simple systems. Moreover, interactions be- 
tween ultracold atoms, characterized by the s-wave scat- 
tering length a, can be tuned via magnetic Feshbach reso- 
nances, giving access to the regime of strong interactions. 
In this letter, we investigate the condensation energy 
of a dilute spin-1/2 strongly interacting Fermi gas with 
variable interaction strength. We show that the conden- 
sation energy can be measured by applying a chemical 
potential imbalance between the two spin states which 
is the analogue of a magnetic field in superconductors. 
We compare our experimental results to a diagrammatic 
theory, finding excellent agreement. 

The experimental setup was presented in [T3]. Our 
system is a quantum gas of 6 Li prepared in a mixture of 
its two lowest energy spin states. The gas is loaded in a 
single beam dipole trap, providing a radial (strong) con- 
finement, while the axial (weak) confinement (z-axis) is 
provided by magnetic coils. This results in a cigar-shaped 
trap. The interactions are tuned using a pair of coils 
in Helmholtz configuration in order to create the large 
homogeneous bias field to tune the scattering length a 
via the 832.18 G Feshbach resonance [TO]. The mixture 
is cooled to quantum degeneracy by lowering the trap 
depth and absorption images perpendicular to the weak 
direction are recorded to obtain the in-situ density distri- 
butions along the z-axis. Previous theoretical [HUE] an d 
experimental studies [T3J [TT] have demonstrated that the 
density profiles of a trapped spin-imbalanced Fermi gas 
can be used to extract the Equation of State (EoS) of the 
corresponding homogeneous system via the pressure for- 

2 

mula: P(/ii, \ii, T) — ^^(ni(z)+fi2(z)), where oj r is the 
radial trapping frequency, and fii(z) — f d r rii(r,z) is 
the doubly-integrated density distribution of spin-species 
»(* = 1,2). 

At unitarity, where the scattering length a diverges, 
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FIG. 1: (Color online) Reduced pressure h(t,b) = 

P(fJ,i, A*2, T)/2Po(p.) of the spin-1/2 unitary Fermi gas, where 
Pq is the T — Fermi pressure of an ideal gas, t the reduced 
temperature ksT/ [i, and 6 = ie unpolarized gas. Black 
points are from from |TJ[]. The Fermi liquid fit is shown as 
a continuous red line, and the extrapolated zero tempera- 
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ture pressure of the normal state £ N is signaled by a red 
cross. Data from [15] is shown in empty blue squares, the 
corresponding fit in dashed blue line, and the extrapolation 
at t = in empty blue square. Bold diagrammatic Monte 
Carlo calculation |16| is shown in solid green line, the ver- 
tical dashed line signals the normal/superfluid transition at 
t c = 0.33. 



we previously measured the pressure of the spin-balanced 
gas as a function of the reduced temperature t = ksT/ii 
(where 2[i = fix + /x 2 ) [T3"] . as well as the pressure of 
the spin-imbalanced gas at t k as a function of the 
spin-polarizing field 6 = ^~^ 2 . We suggested that the 
low-temperature properties of the normal phase of the 
Fermi gas were consistent with a Fermi-liquid behavior 
[17j . As a result, the low-temperature and low-imbalance 
limit of the pressure of the unitary gas can be written as: 
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where -Po(m) = 15^2 {w) ^^ ^ s tne i^eal F erlm S as 
pressure. The response coefficient to temperature t is 
the dimensionless specific heat cy, while the response to 
the polarizing field b is the dimensionless magnetic sus- 
ceptibility x (respectively equal to 5tt 2 /8 and 15/4 for 
the ideal Fermi gas). The magnetic susceptibility has 
been the subject of a previous work [17] . and we focus 
here on the measurement of pressure of the normal phase 
£ N in the t — and 6 = limits. In the (t, b) plane, 
our measurements of the EoS of the unitary gas have 
been performed along two directions: the unpolarized 
gas as a function of temperature h(t, b — 0) (FigjTj), and 
the low-temperature polarized gas versus chemical po- 
tential imbalance h(t — 0,6) (Fig|2ja)). The quadratic 
behavior of the pressure versus both b and t supports 
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FIG. 2: Pressure of the spin-imbalanced gas in the BEC- 
BCS crossover at i = 0. The position of the first-order phase 
transition to the superfluid is shown by a vertical dashed black 
line, (a) Unitary limit. The Fermi liquid fit is shown as a solid 
red line, while the t — Equation of State in the superfluid 
phase is in blue solid line. The pressure of the non-interacting 
gas is displayed as a dotted black line. The t = and 6 = 
extrapolation of the normal phase pressure is shown as a red 
cross, and the condensation pressure as the double-arrows, 
(a), (b), and (c): results of the ladder approximation for the 
normal phase are shown in green for 8 = 0, —0.58 and +0.2 
respectively. 



the Fermi-liquid interpretation of the low-temperature 
thermodynamic properties of the normal phase. How- 
ever, the system will ultimately undergo a second order 
phase transition to a superfluid state and, below a tem- 
perature t c ~ 0.33, the pressure of the spin-balanced gas 
deviates from the t 2 behavior. In contrast at t — 0, the 



spin-imbalanced gas (/Zi 7^ /12) undergoes a first-order 
phase transition to an unpolarized superfluid phase when 
/is (0,0) = /ijv(0, 0) + xb 2 /2. This condition is the ana- 
logue of Eq(2] and at unitarity it yields the critical chem- 
ical potential imbalance b c sw a/0.8, see Fig[2[a). This 
is testified by the discontinuity in the slope of h vs b 2 . 
From Eq.©, and extrapolating the Fermi- liquid behav- 
ior to the zero-temperature and spin-balanced limits, we 
measure the T = dimensionless pressure of the spin- 
balanced unitary gas in the normal phase £^ ' . In the 
first limit (t ->■ 0, b = 0) we find £ N = 0.51(2), while in 
the second one (t = 0, b — > 0), we extract £jv = 0.53(2) 
(see the red crosses in FigfT] and Fig|2|a)). The agree- 
ment between these values, taken for two very different 
limiting regimes, is remarkable and further supports an 
accurate description of the normal phase of the unitary 
gas as a Fermi liquid. This value is in good agreement 
with values calculated using Monte Carlo methods 0.54 
[15] . 0.56 Pi], and 0.52 20 . It is somewhat larger than 
a recent experimental determination, 0.46(1) [15] . 

The problem of the zero temperature balanced super- 
fluid Fermi gas has been the subject of thorough theoreti- 
cal investigations [21]. However, much less work has been 
devoted to the equation of state of the zero temperature 
normal phase [IS]. We show below that our experimen- 
tal results can be quantitatively reproduced using the 
ladder approximation [22j [23] . This theory includes the 
repeated two-body scattering between particles 1 and 2 
described by the scattering length a. In particular, for 
a -1 > 0, it contains the physics of a molecular state. We 
use the finite temperature formalism and take the zero 
temperature limit. The self-energy for the particles 2, 
which physically describes the effect of interaction be- 
tween particles is given by (we take h = 1) 
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where the two-particle vertex T is given by 



r(K,ny 



ni 

An a 



FL(K,£l) 



(5) 



where H(K,fl) is the pair bubble [23]. At zero tempera- 
ture, H(K, fl) can be calculated analytically. The pairing 
instability, signaling a second order phase transition, is 
found using the Thouless criterion r _1 (0,0) = 0. For- 
given /xi and a, this happens for a critical value of the 
chemical potential /i2 C of particles 2. In order to stay 
in the normal phase, we have performed our calculations 
for /12 < fJ-2c- The integration on Q can be performed 
by deforming the integration contour in the half-plane 
5ft($l) < 0. In this way, we pick the singularities of the 
integrand in EqHland get three contributions correspond- 
ing to the pole of (fi — iuj + fii — 2 7„ )~ 1 (^l), the 
branch cut of T(K, 51) (£r) and the molecular pole £Iq(K) 



(for or 1 > 0) of T(K,Q) (E m ) [53]. Sl (K) + 2^ repre- 
sents physically the energy of a molecule of momentum 
K. We find that in the normal phase ilo(K) > 0. As 
a consequence, when we deform the integration contour 
in 3?(J1) < 0, we do not get any contribution from the 
molecular pole of T, and therefore we have S m = 0. This 
is consistent with the physical argument in favor of the 
absence of molecule in the normal phase. Indeed, if we 
had some molecules in the system, they would be con- 
densed at zero temperature. Therefore the system would 
be superfluid, and we are no longer entitled to use Eq|4] 
We deduce the minority density n 2 using the Fermi liquid 
type relation due to Landau 
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where by definition kp 2 = (67r 2 U2) 1 ^ 3 , is the Fermi wave 
vector of particles of type 2. For given /Zi, H2 and a, 
this is an implicit equation for kp2, hence n 2 . As it 
was found in Ref . [24j . we find a non zero density 712 for 
a chemical potential [ii larger than the polaron !24 [25 
chemical potential /i p (/ii). In practice, we fix [i\ > 0, 
then we solve Eqj6jfor a given ^2 > ^ p (/zi). The pressure 
is determined by integrating the density using the Gibbs- 
Duhem relation 

P(m,n 2 ) = Po(m) + f£ d^ [k F , 2 (nun' 2 )} 3 (7) 

For a fixed /ii, we calculate the minority density for in- 
creasing minority chemical potential between /x p (/ii) and 
fi2- For a sufficiently large chemical potential difference, 
the system is normal (the pairing susceptibility does not 
diverge). For sufficiently low b, we calculate the dimen- 
sionless equation of state h(S,b), where 8 is the grand- 
canonical interaction strength, S = h/y / 2m(j,a. For all 
values of 8 < 0, we find a linear behavior of h as a 
function of b 2 . The comparison between experiment and 
theory is shown for 8 = (Unitary limit), 8 = —0.58 
(BCS side of the crossover) and 8 — 0.2 (BEC side) in 
FigsJ2JJa)J2^b) and [2^c) respectively. The agreement is 
very good. However, for increasing a -1 > 0, the values of 
6 in the normal phase become larger and larger, and as a 
consequence the linear fit of h as a function of b 2 , valid at 
low 6, is worse. Still, for 8 = 0.2 the experimental equa- 
tion of state h(S, b) is in good agreement with the ladder 
approximation calculation above 6 C , curved green line in 
FigMc). Within the ladder approximation we have de- 
termined the critical spin polarizing field b c at which a 
pole appears in the vertex function T at zero frequency 
and zero wave vector (Thouless criterion). We found that 
b c was always smaller than the experimental value of the 
first order transition. Our calculation is therefore free 
of any instability singularity in the normal phase. For 
the spin-susceptibility, we also find a good agreement be- 
tween the ladder approximation, experiments, and Monte 
Carlo simulations of IT71. 



Gathering the results from Fig. 2, we now extract the 
zero-temperature dimensionless pressure hjq of the nor- 
mal phase as a function of 5 [17]. The resulting EoS of 
the normal phase h^(8) is plotted in Figp5]as red empty 
squares together with the ladder approximation calcu- 
lation (green line), showing excellent agreement in the 
explored crossover. For comparison, the previously mea- 
sured equation of state of the low-temperature gas in the 
superfluid phase h$ (5) is shown as blue points and a blue 
line fit [T3]. The difference between the superfluid and 
normal pressure at T = thus represents the condensa- 
tion pressure. The superfluid pressure is higher than the 
normal phase pressure, hs(S) > /ijv((5), hence the grand- 
potential is lower and the superfluid state is the stable 
phase at low-temperature. Turning to the canonical en- 
semble the superfluid and normal phase energies £s and 
£,n as a function of the canonical interaction strength 
1/kpa can be computed from the pressure measurement 
of Fig pausing a Legendre transform [25] . The measured 
condensation energy £jy — £s is shown as a solid black 
line in Figj4] 

In the BCS regime, the condensation energy E c can 
be explicitly calculated from the energy of the super- 
conducting and normal states, yielding the well known 

Q A 2 

result E c — %N^j-, where A is the single-particle ex- 
citation gap, and Ep the Fermi energy. Since E = 
lNE F S, a {l/k F a) (where a = S,N), the BCS equation 
becomes: 
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FIG. 3: (Color online) Pressure of the normal /ijv, red empty 
squares, and superfluid hs, blue circles, phases at low temper- 
ature in the BEC-BCS crossover measured in [14] . The green 
line is the result of the ladder approximation. The solid blue 
line is a guide for the eye, while the red solid line is the result 
of Fixed-Node Monte Carlo calculations 17 . The difference 
between the blue and red/green lines is the condensation pres- 
sure. 



Strictly speaking, this formula is valid only in the 
weakly attractive limit A —¥ 0. For an arbitrary interac- 
tion, the condensation energy is given by a more involved 
function of the gap and, based on dimensional arguments, 
it should be written as 
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where F is a (yet) unknown function with F(0) = 1 to 
satisfy the BCS prediction. In the spirit of Landau's the- 
ory, the U{1) invariance suggests that F can be expanded 
with (A/E F ) 2 and as such, the first beyond-BCS correc- 
tion should be proportional to \A/Ep\ 2 . At unitarity 
where A ~ 0.5Ep [TH] this leads to a moderate 25% cor- 
rection to the BCS prediction, which suggests that the 
range of validity of Eq. [8] should extend beyond the strict 
weakly interacting limit [20] . 
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FIG. 4: (Color online) Relation between the condensation en- 
ergy and the superfluid pairing gap. Main panel: measured 
dimensionless condensation energy £n — fs versus interaction 
strength 1/kpa (solid black line). In red points, BCS expres- 
sion (|8| using the values of A measured in [27]. For com- 
parison the prediction from mean-field BCS theory is shown 
in dashed blue line. A Fixed node Monte Carlo calculation 
[17] coincides with the solid black line. Inset: ratio of the 

condensation energy £at — £s to | ( jjr- ) . 

In order to test the BCS expression (l8| in the BEC- 
BCS crossover, we compare our measurement of the con- 
densation energy to | ( ■£- ) using the values of A mea- 
sured by radio-frequency spectroscopy in [27], red points 
in Fig]4] The agreement seen in Fig[4]indicates that, even 
in the strongly interacting regime, the BCS expression 
is remarkably valid. A more stringent test is provided 
by plotting the ratio between left-hand and right-hand 
side of Eq.([8]) (inset of FigQ and we indeed find a ratio 
close to unity. Note that calculating this ratio using BCS 



mean-field theory provides a reasonable estimate (dashed 
blue line in inset of Fig. El) even though the absolute val- 
ues of the condensation energy (dashed blue line in the 
main panel of FigR}, or of the pairing gap are both quan- 
titatively inaccurate in the strongly interacting regime. 

In summary, we have measured the condensation en- 
ergy of a two-component Fermi gas with tunable interac- 
tions. The temperature and spin-polarizing field depen- 
dence of the normal phase pressure are in good agree- 
ment with a Fermi liquid description. A simple ladder 
approximation calculation quantitatively reproduces ex- 
perimental data at zero temperature in the normal phase. 
Future work will explore the critical region and search for 
exotic phases such as the FFLO phase |21) . 
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